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Abstract 

A complete next-to-leading order QCD calculation of deep inelastic electron-photon 
scattering including direct and resolved real photon components is presented. Soft 
and collinear singularities are extracted using the phase space slicing method. Ap- 
plication of the results for the prediction of single- and dijet cross sections at the 
LEP e~'~e~ collider are presented, using the Snowmass jet definition. The depen- 
dence of the cross sections on the transverse momentum and on the rapidities of the 
jets are discussed. 
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1 Introduction 



The production of jets with large transverse momenta pt in 77 reactions with both pho- 
tons being quasi-real has been studied at TRISTAN for quite some time and more 
recently at LEP 0. The presence of a large scale in the process allows perturbative QCD 
calculations. These calculations are available at next-to-leading order (NLO) in QCD 
0, ^, 1^, 1^ and the agreement of the predictions with experimental data is good (see e.g. 

i)- 

In order to reliably understand the hadronic structure of the photon also the process 
should be studied where one of the incoming photons has a large virtuality P, |9|, p!0| . 
The scattering of virtual on real photons, which is the analog to deep inelastic electron- 
proton (ep) scattering, can be achieved in e"'"e~-colliders by a single-tag experiment, i.e., 
one of the leptons is detected under a certain angle and the other escapes unobserved in 
the beampipe. Single-tag experiments have recently been performed at the LEP storage 
rings and the results of these measurements are expected in the near future Some 
studies of deep inelastic electron photon jet cross sections in leading-order (LO), i.e. C(a^), 



have been done in [|T2[. A NLO QCD calculation of jet cross sections in 7*7 collisions 
for relatively small virtualities up to 10 GeV^ has been performed in |[T^, but without 
considering the longitudinal photon polarization. 

The simplest diagram contributing to the deep inelastic electron photon jet cross sec- 
tion in LO is the quark parton model (QPM) diagram, see Figure 1. The real photon 
and the virtual photon couple directly to the charge of the bare quark, producing two 
final state jets. This so-called direct process can be distinguished from resolved processes, 
where the real photon serves as a source of partons (quarks and gluons) that interact 
with the virtual photon. In addition to the high jets from the hard scattering, the 
resolved contribution is accompanied by a "spectator jet" (or photon remnant jet) with 
small Pt- It has been discussed in that both contributions are of the same order in 
the strong coupling constant. Note, that the contribution involving the resolved real pho- 
ton is traditionally called single-resolved in 77 scattering. We will call this contribution 





Figure 1: Examples of LO processes in deep inelastic electron photon scattering from 
e"'"e~ collisions. Left: direct contribution (QPM diagram). Right: quark initiated resolved 
process. 



simply resolved for ease of writing. Similar to ep-scattering, one motivation for studying 
deep-inelastic e7-scattering is to obtain information on the partonic, especially gluonic, 
structure of the real photon from the resolved contribution of the cross section. Fur- 
thermore, the interplay between the hadronic and the point-like part of the real photon 
structure function can be studied. 

As is well known, only qualitative conclusions can be drawn from LO calculations, 
since the results strongly depend on the factorization scale, especially in the resolved 
contribution. Furthermore, no dependence on jet definitions, present in the experimental 
data, can be seen in the LO results. Only the NLO cross sections, which contain an addi- 
tional parton in the final state, depend on the jet definition. The NLO QCD corrections 
to the QPM diagram, i.e., the 0{a'^as) contributions, consist of the radiation of a gluon 
in the final state and of the virtual corrections. The quark-antiquark pair emitted from 
the real photon can produce a collinear singularity in the real corrections, which has to be 
factorized and absorbed into the real photon structure function. In this way, the direct 
and resolved contributions achieve a dependence on the factorization scale My. Thus the 
direct and resolved cross sections are not independent from each other in NLO QCD and 
the clear distinction mentioned above for the LO case does no longer hold. The NLO 
corrections to the subprocesses of the resolved contribution are equal to those obtained 
in deep inelastic ep-scattering, where the photon interacts directly with the partons from 
the proton. The singularities in the resolved contribution that do not cancel after adding 
the real and the virtual corrections are absorbed into the photon structure function, as 
for the direct case. In this paper we will concentrate on the calculation of the 0{a'^as) 
terms in the direct part of the cross section for both, the transverse and the longitudinal 
polarization of the virtual photon. For separating the singular phase-space regions we 
apply the phase-space slicing method [T^. The resolved contributions are taken from 
the ep-scattering case, which has been calculated using the phase-space slicing method in 

MM- 

The paper is set out as follows. In section 2 we discuss cross sections in deep inelastic 
67 scattering in general and review the LO cross sections. We proceed in section 3 by 
performing the NLO calculation of the direct contribution for both polarizations of the 
virtual photon. The techniques for doing the calculations follow essentially those layed 
out in 1^, |18[. Some results for inclusive single- and dijet cross sections are presented in 
section 4 for kinematic conditions and energies encountered at LEPl and LEP2. Section 
5 contains the summary and an outlook. The results of the analytic calculations are 
presented in the appendix. 

2 Leading Order Cross Sections 

2.1 General Structure of the Cross Sections 

To fix the notation we start by writing down the process of jet production in e'^e" scat- 
tering: 

e+{ka) + e-{h) ^ e+{k'^) + e- + Jets + X . (1) 
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We are interested in the case where one lepton radiates a virtual and the other a real 
photon. Of course, it does not matter which of the leptons radiates the virtual photon, 
but for definiteness we suppose this to be the positron. Thus, the subprocess we have to 
consider is 7a(ga) + IbiQb) Jets + X, with qa = ka — k'^, Qb = kb — fc^ and the virtualities 
= —ql and = — = 0. The electron-positron center-of-mass (cm.) energy is 
sh = {ka + k}jY . The hadronic, i.e., 7*7, cm. energy is W"^ = {qa + qbY- Furthermore we 
define the variables 

lakb , qbka , E'^ 

where and E'^ are the energies of the incoming and outgoing electron in the e~^e~ center- 
of-mass system (c.m.s.), respectively. The variable yb gives the momentum fraction of the 
real photon in the electron. 

The cross section da^+e- for the process described above is given by the convolution 



e"'"e" 



J dxbdyb dae+k fk/jixb) F^/e-{yb) ■ (3) 



Here, F^/e-{yb) describes the spectrum of the real photons emitted from the electron 
according to the Weizsacker- Williams approximation |jl9l, which in its' simplest form 
reads 

27r yb V ^1 

The electron mass is rUe and 6max is the maximum scattering angle of the electron. The 
function fk/'y{xb) is the parton distribution function (PDF) of the real photon which 
describes the probability to find a parton with momentum pb = Xbqb inside the real 
photon, where Xb G [0,1]. The direct process is included in formula (|^) by viewing the 
photon as a parton with momentum fraction Xb = 1- It is obtained by summing not only 
over the quark flavors and the gluon, but also over the photon with f-y/.y{xb) = 5(1 — Xh). 
Last, daf,+k gives the positron-parton cross section, which is given by 

where L^^, = i^ka^k'^^, — k'^^kau — g^vkak'^) is the lepton tensor of the positron. In the 
resolved case, H^" is the hadron tensor familiar from deep-inelastic ep-scattering. In the 
direct case, H^^'^ stands for the appropriate tensor describing the photon-photon scattering 



amplitudes (for a definition see e.g. P, ^). Note, that the lepton tensor of the electron 
radiating the real photon integrated over the azimuthal and polar angle of the outgoing 
electron gives the Weizsacker- Williams formula (^. The n-particle phase space of the 
final state particles of the subprocess is rfPS*^"^ and the positron phase space is given by 

^ d<p dygdQ^ 
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The azimuthal angle of the outgoing positron, 0, can be integrated out in the hadronic 
c.m.s. Using the definitions of the trace part of the hadron tensor Hg = —g^^'^H^j.y (denoted 
by g-part in the following), the longitudinal part Hi = {AQ'^) / {s HyaXbYVbVhH ^j.u and the 
unpolarized part Hu = \{Hg + Hl), we find 

1 rr,. _ l + {l-yaf ^ , 4(l-y,) + l + (l-y,)^ ^ 

1 + (1 - VaY ^ , 2(1 - y,) ^ 

yi yi 

The expression Hu is used here to rewrite the result in a form familiar from DIS in ep- 
collisions. It should be noted that the real photon is unpolarized. The longitudinal cross 
section of the real photon is proportional to and vanishes since = 0. 

2.2 Results in Leading Order 

The LO cross sections are (per definition) finite. However, we will use dimensional reg- 
ularization for the NLO corrections and thus we state also the LO results in c? = 4 — 2e 
space-time dimensions. We will neglect all quark masses in the following. 

The Born approximation for the production of two final state quarks in the direct 
process is given by the process 7*7 qq. The resolved gluon and quark initiated subpro- 



cesses 7*g — > qg and 7*(7 qq can be found e.g. in |]T^. Note that the photon induced 
process 7*7 — > qq is very similar to the one where the real photon is replaced by a gluon. 
The photon induced cross section can be obtained from the gluon induced one by keeping 
only the Abelian terms and setting Nc = 1 (hence = 0) for the virtual photon vertex. 
In this way also the NLO corrections can be deduced. 

Using the usual definitions of the Mandelstam variables s = {pb + g^)^, t = {ph — piY 
and u = (pb — ^2)^5 where pi and p2 are the momenta of the final state particles, the 
two-body phase space reads 

,p.(2) _ 1 (Ms + QT] ' 1 dt 

The matrix elements for the QPM process 7*7 —>■ qq read 

H^l*! - qq) = S2aVQtNc Tj^^s, t, u) . (9) 
For the gluon induced resolved matrix elements one finds 

Hfj-f*g-^qq) = 16aayQ^^T^Js,t,u) , (10) 
whereas the quark initiated resolved matrix elements are given by 

Hf,Lh*Q ^ Q9) = 32aasn^Q',CF T^^is, t, u) . (11) 
The definitions of Tgj^^Tgj^ and Tgj^ can be found in Tab. |l|. 
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Table 1: Definition of tlie LO matrix elements j^^T^ ^ and Tj^ in terms of Ti,T2 and 
Ts, defined in tlie appendix. 





photon induced 


gluon induced 


quark induced 


g-part 


T^is,t,u) = Ti{s,t,u) 


T^is,t,u) = Ti{s,t,u) 


T^ = -T,iu,t,s) 


longitudinal part 


T2{s,t,u) = nis,t,u) 


Tl{s,t,u)=Ts{s,t,u) 


Tl = T2{s,t,u) 



3 Next-to-Leading Order Cross Sections 

The NLO corrections are calculated with the help of dimensional regularization. The 
ultraviolet singularities in the one-loop contributions are regularized and subtracted in 
the modified minimal subtraction (MS) scheme. The infrared and collinear singularities 
from the real corrections are likewise calculated in dimensional regularization. 

The 0{as) perturbative QCD corrections to the direct Born process 7*7 qq are 
given by the gluon bremsstrahlung process 7*7 qqg (see Fig. 2, upper row) and by the 
one-loop contributions to this process (see Fig. 2, lower two rows). The NLO corrections 
to the resolved processes constitute the one-loop contributions to the Born terms of section 
2.2 and the four three-body contributions 

7*g qqq , 7*g qq'q' . 

The NLO resolved parton cross sections are of order aa"^, but, after they have been folded 
with the PDF's of the real photon, of the same order as the NLO direct cross sections 



14 1 . In the following subsections we will consider only the direct part in some detail. The 



NLO resolved matrix elements are taken from Graudenz WT\ . 



3.1 Virtual Corrections up to 0{a'^as) 

We deduce the virtual correction to 7*7 qq as mentioned above by keeping the Abelian 
term of the corresponding gluon induced process in order aa^ and adjusting the color 
factors and the coupling. The one-loop process 'y*g — > qq can be calculated by crossing 



from the well-known one-loop corrections in e~^e — > qqg [0, |22|, which has been done in 

Following the above procedure we obtain the results for the g-part and the longitudinal 
part: 



^,^L(7*7-gg) = ^^(^^f^) {4naQfr2NcC,V,AsXn) . (12) 



The expressions Vg and Vl are stated in the appendix. The two-body phase-space is again 
given by eqn (||). 
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Figure 2: The 0{as) corrections to the QPM diagram of Fig. 1. Upper row: real gluon 
emission. Lower two rows: virtual corrections. Additional graphs are obtained hj s —>■ u 
crossing. 



3.2 Real Corrections up to 0{a^as) 

The real corrections from the gluon bremsstrahlung process contain singular regions of 
phase-space that lead to two classes of singularities after integration. The final state gluon 
can become soft or collinear to one of the other final state particles, which will produce a 
final state singularity. Due to the on-shellness of the real initial state photon, the splitting 
term '~f ^ qq will produce an initial state singularity. 

The method used here for handling these regions is the phase-space-slicing method 



15| , where a cut-off ys is introduced in the phase-space integration to separate finite and 
singular regions. The finite integrals are handled numerically which opens the opportunity 
to introduce experimental cuts and jet definitions in a flexible way. The singular integrals 
are done analytically in c? = 4 — 2e dimensions. Both, the analytical and the numerical 
contributions will depend on the slicing parameter. Since ys is a non-physical parameter 
the dependence of these parts on yg should cancel in the sum. We will show this explicitly 
for our results later on. 

After adding the 0{a'^as) terms only single poles proportional to the photon splitting 
function Pq^-y will remain. These are absorbed into the photon structure function. We 
now handle the final and initial state singularities separately. 
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3.2.1 Final State Corrections 



We consider the gluon bremsstrahlung process 



7* (9a) + liQb) q{pi) + q{p2) + g{p^) 



(13) 



and concentrate on the case, where e.g. the invariant r = {P2P3) / {QaQb) vanishes and 
produces a singularity. The case pips can be obtained by simply substituting p2 pi- 
The analytical integration of the three body contributions Hg and Hl for the process 
(p!3|) over rfPS'-'^^ in the singular phase-space region is done by separating the three-particle 
phase-space into the two-particle phase-space (iPS*^^) and a singular phase-space (iPS'^''^. 
In the c.m.s. of the final state particles p2 and ps we define the variable b = |(1 — cosO), 
where 9 is the angle between and p^, and z = (qaPs) / {laQb) ■ Furthermore, defines the 
angle between the two planes defined by the particles {p2,P3} and {pi, g^}. We define the 
variables 



s = s = {qa + qb) 



i= (qb-Pi) 



u = {qb-P2- PzY - 2P2P3 



(14) 



from the three-body momenta. For r ^ the separation of the phase-space is given by 

dps(3) = rfps('^)rfps(2) with i, p, o 



with the function 



Ffl- 



r(2 - 2e) N. 



Gpir) 



1 + 0(r) 



and the normalization constants 
1 

Nb= I dh[h{l - h)]-' 







r^(i-^) 

r(2-2e)' 



>sm 



-2e 



4^7rr(l - 2e) 

r^(i-e) 



(15) 



(16) 



(17) 



The two-body phase-space is once again given by eqn (H). Defining the variable 

Uf = min[— t/(s + Q'^),ys], the limits of integration in ciPS'-^'' are given by r G [0,?/f], 
b e [0,1] and G [0,7r]. 

Expressing the 2 — > 3 matrix elements Hg and Hl by the variables s, t, u, r, b and 
in the limit r — 0, we obtain the approximated matrix elements Hg and H[. These are 
then integrated over the singular region of phase-space: 



^'^ 27rr(l-2e) 



'iTT/i^ 



{AnaQiY 2NcCf Fg^L 



The final results for Fg and Fi are listed in the appendix. They contain the IR collinear 
and soft singularities that cancel against those of the virtual corrections. It is essential 
that the singular terms are proportional to the LO matrix elements and that s, t and u 
reduce to the usual two-body invariants s, t and u for the above limit r 0. 
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3.2.2 Photon Initial State Corrections 



In the direct case the real photon can split into a qq pair that gives rise to a coUinear 
singularity if the partons are emitted parallel. The singularity appears when the invariant 
^1 = ilbPi) / (qaqb) vanishes. We define the new variable 

Zb = 7 r G [r]b, 1 , (19) 

[qaqb) 

that gives the fraction of the momentum that participates in the subprocess after a 
particle has been radiated in the initial state. The variable rjh G [0, 1] is connected to zt, 
through T]i, = XbZb. Furthermore, we define the invariants 

s = {p2+ Psf, i = {zbqb - P2f, u = {zbqb - Psf (20) 

from the three-body momenta. All other definitions concerning the phase-space slicing 
remain as in the previous subsection. In the limit Zi — > the variable s reduces to s. In 
the same limit the three-body phase-space separates according to dPS*-'^^ = dPS'-^-'dPS*-^'', 



where 13 



. .,dz, l-z„ Q-y , , ,„,, 



with 

S - ZbQ^ 



GAz^ 



l-Zi- 



0{zr) . (22) 



s(l - Zb) - ZbQ^_ 

The two-body phase-space is given by equation (|^). The integration over c/PS^*"-' with 
Zi G [0, —u/ (s + Q^)], Zb G [rjb, 1] and G [0, vr] is restricted to the singular region of Zi in 
the range Zi G [0, yj] with ?// = mm[—u/ (s -|- Q^), y^]. 

Expressing the matrix elements Hg and with the variables s,t,u, zi, Zb,b and 
and taking the limit zi 0, one obtains the approximated matrix elements and H^. 
These are integrated according to 



^ ^ Vb 

where Ig and II are stated in the appendix. 

3.3 Finite Next-to-Leading Order Cross Section 

The full NLO cross section da^+^J^ can be written as a sum of the |/s-dependent two- and 
three-body cross sections da2 and da^: 

da^+'^^o = da2{ys) + da,{y,) . (24) 
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The ?/<j-dependence in da^+^J^ cancels up to terms of order i/s In" i/s which have been ne- 
glected in the analytical calculations. Thus, t/g has to be chosen sufficiently small for the 
approximations to be valid. The two-body cross section consists of the final state, virtual, 
initial state and Born contributions da^, da^ , da^ and da^: 

dcT2iys) = da^ivs) + da^ + da'{ys) + da"" . (25) 



When the virtual and final state cross sections are added, the pole terms cancel, as can 
be seen by directly comparing formulae (|D) and (|T) with (^2|) from the appendix. The 



simple pole appearing in the initial state singularities Ig and II is proportional to the 
splitting function Pg«_^(zf,). This function appears in the evolution equation of the PDF 
of the real photon as an inhomogeneous (point-like) term. Therefore, the photon initial 
state singularities can be absorbed into the real photon PDF, according to the procedure 
given in [23, 18]. We define the renormalized PDF fq/e- of a quark q in the electron as 



(26) 



where /^/e- is the LO PDF before the absorption of the collinear singularity. The NLO 
transition function is given by 



M2 



(27) 



with Cqi^f{zi)) = in the MS scheme. In the discussed order, Pg^^{z}y) = 0. The initial 
state contribution da\Yl ~^ jsts) for the photon initial state singularity is calculated 
from the unrenormalized cross section da by 

da\-i*-i ^ jets) = rfcr(7*7 ^ jets) -^j dz,T^^l^{z,, M^^)da^{^*q ^ jets) . (28) 

The cross section da^ contains the LO photon-parton scattering matrix elements given 
in section 2, Hgj^{pf*q — >• qg). The factor (47ryU^/M^)'' in eqn (|27| ) is combined with the 
factor [A-Kfi^/Q'^y in eqn (|23| ) and leads to a factorization scale dependent term of the 
form 



- -Pr 



'47r/i^ 



-p 



(29) 



In this way, the subtracted partonic cross section will depend on the scale M^, as does 
the PDF of the resolved photon. 

After this subtraction procedure, all terms in the two-body cross section da2 are finite 
and the limit e — > can be performed. Eqn (|^) then yields a physically meaningful cross 
section. 
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4 Numerical Results 



We now come to a discussion of numerical results for single- and dijet inclusive jet cross 
sections. Before analyzing these physical cross sections, we discuss some checks of the 
analytical results, one being the dependence of the cross sections on the phase-space 
slicing parameter yg. Furthermore we show the dependence of the cross sections on the 



factorization scale. The multi-dimensional integrals are performed using the VEGAS 
Monte Carlo integration routine. 

At LEP, electrons are scattered on positrons with equal energy. Events are selected 
in a certain angular range of the tag and rejected if there is also an antitag above a 
certain energy and scattering angle. This corresponds to the situation, where one photon 
is quasi-real and the other has a finite, non-vanishing virtuality. 

In this paper we will give predictions at LEPl and LEP2 energies for certain character- 
istic values of the photon virtuality namely = iq, 20, 100 and 200 GeV^ At LEPl, 
the electron and the positron both have the energy E^, = Ep = 45.5 GeV, corresponding 
to a c.m.s. energy of y/s = 91 GeV. The energies of the two incoming leptons at LEP2 are 
Ee = Ep = 91.5 GeV, giving a c.m.s. energy of y/s = 183 GeV. The antitag conditions 
are taken into account by using the Weizsacker- Williams formula (^) with Omax = 25 
milliradians in the case of LEPl and with 6max = 33 milliradians for LEP2. We have 
not imposed cuts on i/a and i/h, i.e., ya,yb G [0, 1]. We will calculate the cross sections in 
the hadronic (7*7) c.m.s. For the theoretical studies here we will refrain from taking into 



account the exact tagging conditions from the LEP experiments |TT[ and will concentrate 
on the above stated energies, although also higher energies will become available at LEP2. 
The jets are clustered at LEP using the Snowmass cone algorithm pS, where hadrons 



i are combined into a single jet, if their distance from the jet center in azimuth-rapidity 
space, given by 



Ri,J=^J{V^-VJ)'+{<P^-<l>J)' , (30) 

is smaller than a cone radius R, i.e., for Rij < R. Here, t]j and 0j are the rapidity and 
the azimuthal angle of the combined jet respectively, defined as 

Etj = Et^ + Et2 (31) 
EnVi + ET2V2 

Vj = ^ , (32) 

ET,(j)i + ET2<h ^22^^ 



E' 



Tt 



The radius R is chosen as in the experimental analysis, namely R = 1. 

For the parton densities of the photon we use the NLO parametrization of Gliick, 



Reya and Vogt |2^, transformed from the DIS^ to the MS scheme. The corresponding 
A value is A*^^) = 200 MeV. Since the energies are quite large, we will choose Np = 5 
flavors for our computations. The renormalization and factorization scales are set equal 
to Q. In particular for the larger values of this is a slightly better choice than Et, 
since the minimum transverse energy of Et,^-^^ = 3 GeV is rather small. We could have 
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chosen Et as a scale, especially for the E'r-spectra, but we prefered to take one scale 
consistently throughout the work. The choice of scale only has a marginal effect on the 
results presented in the next sections. 



4.1 Check of the Analytical Results 

From section 3.3 it is already clear that the poles encountered in the real corrections 
cancel against those from the virtual corrections, which is a test of the singular terms in 
the NLO corrections. The remaining photon initial state singularities have been absorbed 
into the structure function of the real photon. As a further test we have checked that 
the analytic formulae given in the appendix reduce to the expressions encountered in 77 
scattering (photoproduction limit, i.e., 0) stated in [Q. 

A test of the i/s dependent terms is the calculation of the two-body and three-body 
cross sections as a function of yg. The dependence should cancel in the sum of the two 
contributions. In the analytical calculation, terms of order i/s In" i/s have been omitted and 
thus the slicing parameter, which is normalized with the hadronic c.m.s. energy, has to be 
chosen much smaller than one. In Fig. 3 a and b we show the single inclusive direct jet cross 
section integrated over Et > 3 GeV as a function of i/s, integrated over the two different 
Q2_regions G [10, 20] GeV^ and G [100, 200] GeV^ for the above stated LEPl 
conditions. The analytic two-body contributions (dashed) are negative and the absolute 
value is quite large for the smaller values of i/s- The three-body contributions (dotted) are 
large and positive. The sum of these (full) is independent of the slicing parameter in the 
shown region of ys, as expected. Although it is theoretically safer to choose smaller values 
of ys, because of the omission of the (9(|/s)-terms, it is from the practical (numerical) 
point of view convenient to have larger values of ys, since the compensation between the 
two-body and the three-body parts are then smaller and the statistical errors are likewise 
smaller. We will choose ys = 5 ■ 10^^ in the further calculations. 

In Fig. 3 we also compare the NLO cross section with the LO result (the dash-dotted 
line), which does not depend on ys. For the LO cross sections we have chosen the same two- 
loop as formula and the same PDF as in the NLO cross sections, so differences between 
the LO and NLO results are due to differences in the LO and NLO matrix elements. The 
NLO curves are smaller than the LO order ones and the corrections are around 10% for 
the smaller region. For the larger region the NLO corrections reduce the Born 
result by approximately 30%. 

The resolved contributions are also independent of ys, which we explicitly checked but 
do not show here. Similar tests of the resolved contribution have already been performed 
for the case where the real photon is substituted by a proton in [O . 



4.2 Scale-Dependences 

Having checked the pole terms, the limiting behavior and the ?/s-independence, which 
shows that the program for the numerical calculations is reliable, we proceed by discussing 
the dependences of the cross sections on the renormalization- and factorization-scales in 
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LO and NLO. This has been done before in the hmiting case of 0, i.e., for 77- 

scattering, in 0. 

In Fig. 4 a and b we show the variation of the NLO direct and LO resolved cross 
sections with the factorization scale (normalized with Q) for the same kinematic 
conditions and regions of virtuality as in Fig. 3. In both regions, the LO resolved cross 
section (dotted line) grows by about 50%, going from = to = 2Q. This 
dependence is due to the factorization scale dependence of the photon structure function 
alone, since the LO partonic cross sections are finite and no initial state singularities have 
to be treated. The dependence of the LO resolved contribution is completely canceled 
by the NLO direct contribution, shown as the dashed line. The sum (full line) shows 
no dependence on the factorization scale in the plotted region. The dependence of the 
NLO direct contribution on the factorization scales comes in through the absorption of 
the singular terms in the initial state into the photon structure function, as explained in 
section 3.3. The resolved NLO corrections to the partonic cross section do depend on 
the factorization scale, in contrast to the resolved LO partonic cross sections. One sees 
however only a small overall dependence on the factorization scale of under 4% by going 
from = |Q to = 2Q, when the NLO direct and resolved contributions are added 
to obtain the full NLO cross section (we do not show these plots here). 

The renormalization scale dependence of the resolved cross section is reduced by going 



from LO to NLO, as has been shown in [131. For the direct cross section, the LO result 



is independent of the renormalization scale, because the strong coupling constant as does 
not appear in the QPM diagram. Thus, the NLO cross section is actually LO in the 
strong coupling. As we will see later on, the NLO correction to the LO direct process 
is small and thus the renormalization scale dependence, which is only seen in the NLO 
corrections, will likewise be small. In addition, the running of the strong coupling is 
compensated by including the virtual graphs in the NLO direct process. We found only 
a small variation of the NLO direct cross section of about 10% between the smallest and 
largest scale considered here. Thus, summing the direct and resolved contributions we 
find a reduced renormalization scale dependence due to the reduced dependence of the 
resolved components. 

4.3 Inclusive Single-jet Cross Sections 

We start this section by discussing the single-jet inclusive r^-distributions 

dEr J,^,, (34) 



drjdQ"^ J dE'rdQ'^dri 

in the region |?7| < 3, integrated over Et > 3 GeV for LEPl conditions. We show the 
direct and resolved cross sections and their sum in Fig. 5 and look at the longitudinal and 
transverse contributions to the individual direct and resolved components in Fig. 6 and 
Fig. 7. 

Fig. 5 a,b,c and d show the direct (dashed) and resolved (dotted) contributions as 
well as their sum (full line) for the four values defined above, = 10, 20, 100 and 
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200 GeV^. The direct contributions are peaked at 77 = and are quite symmetric, which 
is to be expected, since we are in the 7*7 c.m.s. A small asymmetry arises from the 
mass of the virtual photon. The resolved cross section is peaked at rj = 0.5 and shows 
a stronger asymmetry. This is due to the additional integration over the momentum 
fraction of the parton in the resolved photon. The peak is in the positive ?7-region, since 
the virtual photon is chosen to travel along the positive z-axis. Four all values, the 
direct contribution is the dominant one, although the resolved becomes more important 
for larger values and in general for positive ?7's. This shows that the direct contribution 
falls of slightly stronger with rising compared to the resolved contribution. The reason 
for this is the relatively low cut on Et, where the resolved component is expected to be 
more important. 

In Fig. 6 a,b,c and d we show the direct cross section, split up into the transverse and 
longitudinal parts, compared to the full LO cross section. The transverse part (dashed 
line) is denoted by and is obtained from the full cross section by setting Hl = 0, so 
only the g-part is in ar- Likewise, the longitudinal part ctl (dotted line) is obtained by 
setting Hg — 0. The sum of both gives the full direct contribution, denoted here as an (full 
line). For the plots a and b, the transverse part is dominant, although the longitudinal 
part is important, especially in the central region. One observes a plateau around 77 = 
for the transverse part, which becomes a dip for the two larger Q^-values in Fig. 6 c and 
d. This local minimum is compensated in the full direct cross section ctd through the 
longitudinal component which has a maximum around 77 = 0. From this one sees that 
the longitudinal polarization of the virtual photon gives an important contribution to the 
direct cross section. The NLO direct cross section is compared to the LO cross section 
(dash-dotted line), containing both polarizations of the virtual photon. In Fig. 6 a and b, 
the LO and NLO cross sections can not be distinguished in the logarithmic plot, but the 
NLO curves lie around 5-10% lower than the LO ones, which is consistent with the results 
shown Fig. 3. For the two larger Q^-values, the NLO curve lies approximately 25% below 
the LO cross section in the central ?7-region. Furthermore, the NLO cross section has a 
slightly stronger fall-off towards the edges of the 1] phase space than the LO cross section. 

The dip observed in the NLO transverse cross section is already present in LO, which 
we do not show here. It is not present in photoproduction and for the smaller values. 
In LO the reason for the local minimum at larger is twofold. First, the interference 
term in the Born matrix elements —Q^s/{ut) is negative and sharpley centered around 
77 = 0. For small this contribution is negligable, but going to larger virtualities the 
subtraction of the interference term in the central 7;-region leads to the observed local 
minimum. Second, the individual terms t/u and u/t va. the Born matrix elements are 
peaked on opposite sides, shghtly away from 77 = 0. For = their sum still gives a 
curve with an absolute maximum at 7/ = 0. At larger values though the maxima of the 
terms t/u and u/t move to larger (smaller) 7/- values, so that their sum leads to a slight 
local minimum near 77 = 0. These effects, observed in LO, have their correspondences in 
NLO. 

Next, we consider the resolved components, again separated into transverse and lon- 
gitudinal parts, in Fig. 7 a,b,c and d. The curves are labeled as in Fig. 6. The transverse 
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part of the cross section is dominant for all values and the longitudinal part even 
becomes less important for the larger virtualitics and is one order of magnitude smaller 
compared to the transverse part at = 200 GeV^. The longitudinal part is peaked 
more in the positive rj region for the smaller values than the transverse part and it 
is shifted to the smaller jy's for the larger Q'^ values. Here, also the LO curves and the 
NLO curves are not distinguishable in the logarithmic plots of Fig. 7 a and b. For the 
resolved cross section, the NLO cross sections are around 10% larger than the LO cross 
sections, in contrast to the direct case, where the NLO cross sections lie lower than the 
LO cross sections. It should be mentioned that the transverse part of the gluon induced 
resolved cross section also shows a dip near 77 = 0, for the same reasons as in the case 
of the transverse direct cross section. However, the gluon induced process gives only a 
rather small contribution to the resolved cross section and the transverse part of the quark 
induced cross section does not have a local minimum in the ry-distribution, since here the 
interference term is positive. 

In connection with the £'r-distributions it will turn out to be important that the direct 
cross sections have a larger plateau in the central ry-region than the resolved components, 
which are peaked more sharply aroTind rj = 0.5 and fall off stronger towards the edges 
of the 7] phase space. The effect is stronger for the larger values. This is due to the 
special form of the transverse direct cross section, which is relatively broad already for 
the smaller values and has two maxima at the edges of the rj spectrum for the larger 

values. We further mention that the rj distributions become narrower when integrated 
over Et > Et for larger values of Er due to the stronger kinematical restrictions in 
the region of larger Et- This holds for both, the direct and the resolved components. 

We now come to the E^ distributions of the single-jet inclusive cross section 



for the same Q^- values as in Fig. 5. The rapidity is integrated over the central region 
\ri\ < 2 in the hadronic c.m.s. We show the NLO distributions of the direct (dashed) and 
the resolved (dotted) component of the cross section and the sum (full line) of the two in 
the E^-range Et G [3, 11] GeV in Fig. 8 a,b,c and d. 

As can be seen, the direct component is the dominant one in the whole E'^-range for 
the two smaller Q^-values. In addition the resolved contribution falls off stronger with 
rising Et than the direct component. The stronger fall-off of the resolved component 
holds for the two larger Q^-values as well, but here the resolved component becomes 
comparable to the direct contribution for the smallest Et values. The flattening of the 
Et distribution of the direct component for the small -Br's, already visible at = 100 
GeV^, even leads to a crossing of the direct and resolved cross sections at £'t = 3 GeV 
for the largest Q^-value, Fig. 4d, so that the resolved is larger than the direct. 

This behavior is a purely kinematical effect and stems from the fact that the 77- 
integration is restricted to \rj\ < 2. We have mentioned above that the direct cross 
sections have a broader r^-spectrum than the resolved ones. Cutting on r] thus leads to 
a stronger cut on the direct than on the resolved components. The effect becomes more 




(35) 
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important for the larger Cf" values because the direct contributions have an even broader 
?7-spectrum for = 100 and 200 GeV^ than for the smaller values, as already dis- 
cussed above. Since the r^-spcctra become narrower for the larger -Et's, the cut on 77 does 
no longer have such a large effect on the direct cross section and the direct component 
again dominates over the resolved component for the larger E't's. 

As a main result we find that the resolved contribution is small compared to the di- 
rect contribution and therefore there is little hope to learn about the parton distributions, 
especially the gluon distributions, in the photon from deep-inelastic e7-scattering. Fur- 
thermore, the resolved component is largely suppressed for larger transverse energies E't, 
which agrees with the expectation that the point-like part of the real photon is dominant 
at large scales. 

4.4 Inclusive Dijet Cross Sections 

We now come to the presentation of inclusive dijet cross sections for LEPl and LEP2 
energies. In Fig. 9 a,b,c and d we show the dijet cross section da"^^^^ / dET^dQ"^ integrated 
over r/i, ?72 G [—2, 2] as a function of the transverse momentum of the trigger jet Et^ for 
the same Q^-values as in Fig. 5 (a)-(d) for LEPl conditions. We see much the same 
behavior as in the case of the single-jet cross sections, namely that the direct component 
is dominant and the resolved one falls off stronger with rising Ex- As in the single-jet 
case, the direct cross section has a flattening towards the smaller E^s for the largest Q"^ 
values, which seems to be even stronger in the two-jet case. The reasons for the flattening 
are the same as discussed above in the single-jet case, only here we have a cut on the ?7's 
of both jets, which increases the effect of the 77-cuts and thus the flattening of the Et 
distributions. 

In Fig. 10 a,b,c and d the dijet cross section integrated over 771,772 G [~2, 2] is shown 
in the range Et G [3, 11] GeV for LEP2 conditions, i.e., for slightly different antitagging 
conditions than for LEPl and for larger energies. Therefore, the cross sections are larger 
than for the dijet cross sections for LEPl. All other conclusions remain unchanged. The 
phenomenon of the flattening turns out to be quite obvious here. As can be seen in 
Fig. 10 d, the direct contribution even has a slight maximum in the Et spectrum for 
Et — 4: GeV and by going to smaller Et values the direct cross section becomes smaller. 
This is compensated by the rise in the resolved component, so that the sum of the direct 
and the resolved component rises by going to smaller Ets. 

5 Summary and Outlook 

We have presented the calculation of the direct component of jet production in deep in- 
elastic electron-photon scattering in NLO QCD. Transverse and longitudinal polarizations 
of the virtual photon have been taken into account. The singular regions of phase space 
have been extracted with the help of the phase-space slicing method. The initial state 
singularities of the real photon have been absorbed into the PDF of the real photon. Cross 
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sections have been obtained by adding the direct and resolved real photon parts. The 
spectrum of the real photon has been approximated by the Weisacker- Williams formula. 
We have shown that the scale dependence of the NLO cross sections are reduced. We are 
thus in the position to compare NLO calculations to experimental data for deep inelastic 
67 scattering, taken e.g. at LEPl and LEP2. 

We have presented rj- and i^y-distributions of inclusive one- and two-jet cross sections 
for LEPl and LEP2 energies for photon virtualities between 10 and 200 GeV^. The direct 
component gives the dominant contribution to the cross section. The NLO direct cross 
sections are about 20% smaller than the LO cross sections for > 100 GeV^ and about 
the same size for around 20 GeV^. The longitudinal part of the cross section plays an 
important role and can contribute up to 50% in the full cross section. 

It is of course possible to study the influence of different parametrizations of the real 
photon PDF on the resolved cross section. We have refrained from showing such compar- 
isons here, because the resolved plays only a minor role in the cross sections. Furthermore 
it is possible to study the effects of cuts on other observables, such 9jS X .y ^ COS 9* or the 
invariant jet mass, in NLO, as has been done for 77-scattering in |^. However, data on 
these kind of variables are not expected in the near future. 

It has been discussed recently in in connection with ep-scattering at HERA that 
contributions from the resolved virtual photon might contribute to jet cross sections up 
to virtualities of 100 GeV^. The resolved virtual photon contributions have to be taken 
into account when terms of the type InQ'^/E^ become large and have to be absorbed 
into a virtual photon structure function, as is the case for ^ E^. The subtraction 
procedure has been worked out in [^] for the case of ep-collisions and is very similar for 
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In this paper we have not taken into account the effects of a resolved virtual photon. 
However, the case of ep-scattering corresponds to the small resolved contribution in the 
7*7 case. The resolved virtual photon effects in the direct cross section will presumably 
be small for larger values and it is not clear how the double-resolved contribution, 
i.e., the contribution where the real and the virtual photons are both resolved, behaves. 
We will extend our studies presented here in the future to include the effects of resolved 
virtual photons. 
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Appendix 



Born Terms and Virtual Corrections 



The Mandelstam variables s,t and u, as well as the invariants Q'^,SH,Xb and r/a are 
defined as in section 2 in the main part. The LO matrix elements in c? = 4 — 2e space-time 
dimensions are given by 



Ti{s,t,u) 
T2{s,t,u) 

Ts{s,t,u) 



t u\ 2Q'^s 



ut 



2e , 



4g^ 



4Q2 u 



{sHXbyaY 2 



(36) 
(37) 
(38) 



The one-loop contributions to the process 7*7 qq depend on the two-body variables 



s, t and u. The function L{x, y) appearing in the virtual corrections is defined as [jT7 



L{x,y) 



In 



X 



— lim Re 



In 
C2 



y 




In 


X 


Q' 






f X 








+ 11] j 



In 



1-^ 

' y 



In 



In 



I + ^2 ( + 



+ 



6 



(39) 



where C2[x) is the Dilogarithm function. For the g-part of the virtual corrections we find 

27r2 



Ti{s,t,u) 



+ 4 In 



e V2 
2s 

+ 



ut 



, -M / 4s 2m 



V S + t (s + t)2 

s2 + (s + t)2 



-t /4s + 2t 

In — r 



ut 



2L -s,-M 



■ut 



2L(-s,-t) 



\ s + u (s + u)^^ 
s^ + (s + m)2 



+ 



4s s s 
+ + 



,u + t u + s s + t 
For the longitudinal part we find 

'2 2/3 



s s u t 

-+-+-+- 

u t t u 



(40) 



VL{s,t,u) = T3(s,t,u) 



+ T^lllT^ 



-u f 2t 



+ 



e \2 

ut 



-In 



1 -t ( 2u ut 
+ - In — — + 



2 \s + t {s + ty) 2 Q'^\s + u {s + u 



+ 14 + 



+ 



s + u s + t 



L{—u, — s) — L{—t, —s] 



(41) 



18 



Final State Corrections 

The final state corrections to the process 7*7 —>■ qqg depend on the invariant mass cut-off 
i/F and on s, t and u. Neglecting terms of order e, we find 



The longitudinal correction in the final state, Fl, is simply obtained by replacing the 
Born term Ti by T3. 

Photon Initial State Corrections 

The photon initial state singularity depends on the cut-off parameter yj, the additional 
variable of integration zi, and on s, t and u. Again, terms of order e have been neglected: 



The longitudinal term, 7^, can be obtained by replacing the Ti by T3, as for the final 
state corrections, keeping the invariants t and s exchanged. The Altarelli-Parisi splitting 
function in ( ^31) is given by 




(42) 




(43) 




(44) 
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Figure 3: Single-jet inclusive direct cross section da^^^*" integrated over Et > 3 GeV and 
< 2 as a function of?/,, (a) ^ [10,20] GeV^; (b) ^ [100,200] GeV^. 
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Figure 4 a,b: Single-jet inclusive direct cross section da^^^* integrated over > 3 GeV 
and Ir^l < 2 as a function of the normalized factorization scale M^/Q for the same 
intervalls as in Fig. 3. 
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Figure 5: Inclusive single-jet cross section da^^^* / drjdQ'^ integrated over Et > S GeV as a 
function of the rapidity t]. (a) = 10 GeV^; (b) = iq GeV^; (c) = 100 GeV^; 
(d) Q'^ — 200 GeV^. Direct: dashed line; resolved: dotted line; sum: full hne. 
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Figure 8: Inclusive single-jet cross section da^^^* / dExdQ'^ integrated over 77 as a function 
of the transverse momentum Et for the same Q^-values as in Fig. 5 (a)-(d). 
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Figure 9: Inclusive dijet cross section da'^^^* / dEx^dQ'^ integrated over 771,772 £ [~2, 2] as a 
function of the transverse momentum of the trigger jet Et^ for the same Q^-values as in 
Fig. 5 (a)-(d) for LEPl. 
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Figure 10: Inclusive dijet cross section da"^^^* / dET^dQ"^ integrated over 771,772 £ [~2,2] as 
a function of the transverse momentum of the trigger jet Et^ for the same Q^-values as 
in Fig. 5 (a)-(d) for LEP2. 
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